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c/2 | The problem of defining the average kinetic energy of statistical systems is ad- 

dressed. The conditions of applicability for the formula, relating the average kinetic 
energy with the mass derivative of the internal energy, are analysed. It is shown that 
incorrectly using this formula, outside its region of validity, leads to paradoxes. An 
equation is found for a parametric derivative of the average for an arbitrary operator. 
A special attention is paid to the mass derivative of the internal energy, for which a 
general formula is derived, without invoking the adiabatic approximation and taking 
into account the mass dependence of the potential-energy operator. The results are 
illustrated by the case of a low-temperature dilute Bose gas. 

> 

m 

(N 

Key words: Bose-Einstein condensate; dilute Bose gas; kinetic energy; parametric vari- 
! ation of statistical averages 

O 

+-> 

03 

s 

C ' PACS: 03.75.Hh, 03.75.Nt, 05.30.Jp, 05.70. Ce 

O ' 
O 

> 

X 



l 



1 Introduction 



There exists a known formula by Landau and Lifshitz [1] (see §15) connecting the average 
kinetic energy of a statistical system with the mass derivative of the internal energy. This 
formula has recently been used in several papers for evaluating the average kinetic energy of 
low-temperature dilute Bose gas. This gas is known to be well described by the Bogolubov 
theory [2-5]. However, there is an ambiguity in calculating the average kinetic energy in 
the dilute-gas approximation [6]. To avoid this ambiguity, several authors have applied the 
Landau-Lifshitz formula for evaluating the kinetic energy. Unfortunately, one usually forgets 
the conditions, under which this formula was derived, because of which the formula employed 
outside its region of applicability, leads to incorrect results. 

To be concrete, let us recall this formula [1] together with the conditions of its validity. 
One considers a statistical system described by a Hamiltonian 

H = K + W, (1) 

which is a sum of the kinetic-energy operator K and of the potential-energy operator W . 
The kinetic-energy operator is assumed to have the standard form of an operator for non- 
relativistic particles, when K oc l/2m, so that 

am m 

where m is the particle mass. The potential-energy operator is supposed to have no depen- 
dence on the mass, such that 

dm ^ ^ 
The total internal energy of the system is given by the statistical average 

E = <H> (4) 

of the Hamiltonian (1). Let the mass be varied under the fixed temperature T, volume V, 
and entropy S. Then one has the average kinetic energy 



<K> = -m — , (5) 




which is the Landau-Lifshitz formula [1]. Here and in what follows all parametric derivatives 
are taken under fixed temperature and volume, because of which this is not specified each 
time to make the resulting expressions less cumbersome. 

It is, however, important to stress that fixing the entropy implies the adiabatic approxi- 
mation. Really, the entropy, by its definition, is a function S = S(T, m) of temperature and 
mass. Hence, fixing the entropy means that the temperature becomes a function T = T(S, m) 
of the given entropy and mass. Therefore, the derivative (dT/dm)s ^ is nonzero. That 
is, in general, it is impossible to simultaneously fix the temperature and entropy. The latter 
is admissible only for some extremal cases, say, when T = and, hence, S = const. This is 
why formula (5), in the general sense, corresponds to an adiabatic approximation, becoming 
exact solely for those cases, when both T and S can be made constant. 
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But even when T = 0, hence S = const, it is necessary to keep in mind other conditions, 
only under which formula (5) is valid. Going outside the region of applicability of this 
formula can yield principally wrong results. As an illustration, let us consider the uniform 
low-temperature dilute Bose gas described by the Bogolubov theory [2-5] (see also the recent 
review papers [7-9]). The particle interactions are characterized by a scattering length a s . 
We assume that there are no bound states, so that a s does not depend on the mass. At 
zero temperature and in the first order with respect to the interaction intensity, the internal 
energy per particle is the ground-state energy 

N m 

where p = N/V is the average particle density. It is a trivial task to check (see the original 
works [2-5], review papers [7-9], or the textbook [10]) that the average potential energy for 
this case is < W >= E. At the same time, invoking formula (5), one gets < K >= E. 
By definition (4), one has E =< K + W >. Since < W >=< K >= E, one comes to the 
equality E = 2E, implying that 1 = 2, which is, of course, senseless. Such a paradoxical 
conclusion stems from the fact that Eq. (3), which is a necessary condition for the validity 
of formula (5), has been forgotten. 

In the following Section 2, we derive a general equation for the derivative, with respect 
to any parameter, of a statistical average for an arbitrary operator. In this derivation, the 
adiabatic approximation is not involved. The obtained general formula is illustrated by 
considering the derivatives of the internal energy, including the case of the dilute Bose gas. 
In Section 3, we show how the average kinetic energy can be evaluated directly from the 
Bogolubov theory [2-5]. 



2 Parametric variation of statistical averages 

Consider a statistical system characterized by a Hamiltonian H. In equilibrium, the average 
of an operator A is given by the statistical average 

< A > = Tr pA , (6) 

with the Gibbs statistical operator 

-p6 



e 



Tre-f> H 



Throughout the paper, the Heisenberg representation for the operators is kept in mind, and 
the system of units is used, where the Boltzmann and Planck constants are set to unity, 
k B = 1, h = 1. 

Let the Hamiltonian depend on a parameter A, such that H = H(\), and the considered 
operator be also dependent on the same parameter, A = A(X). Taking the parametric 
derivative of two traced operators, say A and B, one should keep in mind the trace property 



Tr AB = Tr BA = - Tr 
2 



A, B 
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where the notation for the anticommutator 

A, b] = Ab + bA 

V J + 

is used. That is, the parametric derivative 



d 1 
^- Tr AB = - Tr 

dX 2 



9 A 



+ 



A 



dB 



has to always be represented in the symmetrized form. 

For what follows, we shall also need the notation for the covariance of two operators, 



cov [A, B) = -< 



A, B 



> - < Ax B > 



This is symmetric, 

cov (i, B) = cov (B, i) , 
and possesses the property of additivity, 

cov (A + B, C)= cov (i, C) + cov (B, C) . 

The diagonal covariance defines the operator dispersion 

A 2 (i) = cov (i, A) = < A 2 > - < A> 2 . 

Taking the parametric derivative of an operator average, we shall always keep the tem- 
perature and volume fixed. But no restrictions on the entropy will be imposed. Then 



9 a dA a (a d ^ 

-<A> = <->-/?covU- 



(7) 



Of particular importance is the case, when the operator A is self-adjoint, representing an 
observable quantity, such that 

(if k ,Aif p ) = (Aip k ,<f p ) . 
Let also the operator A commute with H, thus representing an integral of motion, 



A, H 



= . 



Under these conditions, the operators H and A possess a common set of eigenfunctions, 

Hip k = E k ip k , Aip k = A k (p k . 

Owing to the Hamiltonian parametric dependence H = H(X), the eigenfunctions also depend 
on this parameter, ip k = (p k (\). For the normalized functions, one has 

d 

0^ (<Pk,<Pk) = , (<p k ,<Pk) = l- 



It is straightforward to check that 

9 ( 3 \ ( dA \ 

— [y k ,A(p k j = I <p k , — <p k I . (8) 
Comparing Eqs. (7) and (8), and recalling that 

Hd\ <A> = d\^ A ^ h H < d\ >= [*°> ox^) ' 

with ip being the ground state, we see that for a self-adjoint integral of motion A, one finds 
the limit 

A, w j=o. (9) 

In particular, when A is H, then for the ground-state energy 

E = (v? , Hipo) 

one gets 

dE ( dH \ /rT1 rt . 

and the limit 

(dH\ 
H, -^-1=0. 
dX J 

For the internal energy E =< H > at any temperature T > 0, we obtain the general 
formula 

_ =< _>_ /3cov ^, _) . (10 ) 

Clearly, the second term in the right-hand side of Eq. (10) is not necessarily zero at finite 
temperature T > 0. 

The first term in the right-hand side of Eq. (10) can be expressed through the derivative 

dF /,,, 

ax =< ^x > 

of the free energy 

F — —T In Tr . 

Hence, Eq. (10) can be rewritten as 

dE OF n ( ~ dH\ 

_ = __ />OT ,^,_j. (12) 

Due to the equality F = E — TS, one has 



Then another form of Eq. (10) is 



dE OH ^dS 

9A =< 9A >+T ar (13 > 
From the comparison of Eqs. (12) and (13) it follows that 

The same expression (14) can be derived from the direct differentiation of the definition for 
the entropy S = — TrpTnp, for which we have 

dS ( dH\ 

Substituting here 

lnp = (3(F-H), 
we return to Eq. (14). According to the limit (9), one gets 

hm T ^ = . 

t^o OX 

Let us emphasize that the second terms in the right-hand sides of Eqs. (10), (12), and 
(13), generally, are not zero for finite temperatures. It is only in the adiabatic approximation, 
when one has 

(dE\ (8F\ OH (as n \ 

The situation here is the same as is mentioned in the Introduction. At finite temperatures, 
the entropy is a function S = S(T, A). Hence fixing the entropy requires that the temperature 
be a function T = T(S,X) of the parameter A, so that (dT/d\)s ^ 0. In general, it is 
impossible to fix simultaneously both the temperature and entropy. This could be done 
either in some marginal cases or approximately. 

Let us now take the particle mass m as a parameter A. Then, with the derivative 

5* = -i* + f, (16) 

from Eq. (10) we find 

t— = <K> + < — >+T — . (17) 

om m om om 

This formula is a generalization of the Landau-Lifshitz equation (5). At zero temperature, 
this becomes 

dE 1 ~ dW 

= - - <K> + <^-> (T = 0). (18) 
om m om 
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The second term in the right-hand side of Eq. (18), generally, is not zero. It is only under 
the supposition (3), when Eq. (18) reduces to the Landau-Lifshitz formula (5). 

An obvious case, when the condition (3) is not valid, is exactly the case of the dilute 
Bose gas, whose interaction energy is proportional to the interaction intensity 

$ = 4tt — . 

m 

The potential of pair interactions for the dilute gas is proportional to $o> whether the Fermi 
contact potential or the Huang effective potential are considered [10]. Therefore for the 
dilute gas, one has 

dW 1 - 

^- = -±W. (19) 
dm m 

Note that relation (19) is valid for many systems, but not solely for the dilute gas. Together 
with Eq. (2), this gives 

f = -±H. ,20) 
dm m 

From Eq. (14), we find 

|£ = -f eov U f) = £ A 2 (H) . (21) 
dm \ dm J m 

Remembering the definition of the specific heat 

C ~ 1 ( 8E \ ~ A2{6) (22) 
Cv= N{df) v -l^' (22) 

we get 

^ = - NC V . (23) 
am m 

The latter equation emphasizes again that the adiabatic approximation, in the sense of the 
independence of the entropy from the mass, can become asymptotically exact solely for 
T — > 0, when Cy — > 0. But for finite temperatures, when the specific heat can be rather 
large, one cannot neglect the dependence of S on m, so that dS/dm ^ 0. 
Finally, from Eq. (17) we find 

7T- = - ( NT °v - E ) ■ ( 24 ) 
dm m 



At zero temperature, we have 

where the limit 

is taken into account. 



^ = -- (T = 0), (25) 
dm m 



lim TC V 
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It is formula (24), instead of Eq. (5), which is to be used, when no adiabatic approxima- 
tion is valid and the potential-energy operator depends on the mass according to Eq. (19). 
At zero temperature, Eq. (24) reduces to Eq. (25). For the low-temperature dilute Bose gas, 
for which E/N = 2npa s /m, Eq. (25) becomes a simple identity. In this way, taking the mass 
derivative correctly, without forgetting the dependence on the mass of the potential-energy 
term, results just in an identity, but by no means this defines the kinetic energy. 

The situation diversifies when, in addition to the interaction energy, the potential energy 
term contains an external potential depending on the particle mass. For instance, the exter- 
nal potential can be a confining potential acting on atoms in a trap. Let the potential-energy 
operator consist of two terms, 

W = V + U , (26) 

of which the first term corresponds to particle interactions, with the mass dependence accord- 
ing to Eq. (19), while the second term describes a power-law trapping potential, proportional 
to m n l 2 . For a parabolic trap, n — 2. Thus, we have 

Om m om 2m 

Consequently, instead of Eq. (19), we get 

dW 1 - 1 / n\ - 

JT = - Lw + - 1+ o ) U ■ ( 28 ) 
om m m \ 2 J 

Then relation (17) transforms to 

|£._I £ + I( 1 + |) <£r> + r ^. (29) 

om m m \ 2 J om 
For a harmonic trap at zero temperature, we obtain 

= E + — <U > (T = 0). (30) 

om m m 

This again does not allow one to get in a simple way the average kinetic energy. 

Thus, in any case, considering whether homogeneous systems or trapped atoms, one must 
accurately take into account the mass dependence of the potential-energy term. Then no 
paradoxes appear. While blindly applying Eq. (5) outside of the region of its validity yields 
principally wrong results. 

It is worth emphasizing that Eq. (25) is an identity, for the dilute Bose gas at zero 
temperature, for all known orders of the perturbation theory in powers of the dimensionless 
coupling a = pa z s <C 1. The ground-state energy, as an expansion with respect to this 
coupling, is 

— = 2tt ^— ^ f 1 + bia 1 / 2 + b 2 a + b' 2 a In a) . 
N m v ' 

The coefficients in this expansion, found in Refs. [11-14] and summarized in the review 
article [15], are 

6i = 4.814 , b 2 = 71.337 , b' 2 = 19.654 . 
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As is evident, Eq. (25) is an identity for the ground-state energy of the dilute Bose gas 
in all orders with respect to a. Moreover, the above expansion can be extrapolated, by 
invoking resummation techniques, to arbitrary values of a G [0, oo). Consequently, relation 
(25) seems to be valid, not merely for the dilute Bose gas in the weak-coupling limit, but for 
this gas with arbitrarily strong couplings. 



3 Kinetic energy of Bose-condensed systems 

The direct definition of the kinetic energy per particle, for a homogeneous Bose-Einstein 
condensed system, is 

K=^^ = -y:—n k: (31) 
N N ^ 2m { } 

where n k is the momentum distribution of atoms. The contribution of condensed particles is, 
as is evident, zero because of the zero momenta of particles in the condensate. For a dilute 
Bose gas, the number of noncondensed particles, with a momentum k, is much smaller than 
the number of condensed particles, n k <C No ~ N. The potential energy is mainly due to the 
interaction of the condensed atoms, so that < W >= 2irpa s /m. But the kinetic energy is 
caused by the noncondensed particles, because of which it must be of the higher order with 
respect to the interaction intensity $ — 4:7ia s /m, hence of the higher order in the scattering 
length a s . From here, even with no calculations, it is obvious that the average kinetic energy 
< K > cannot be of the same order in a s as the ground-state energy E W >. Recall 
that incorrectly employing Eq. (5) one would get < K >^ E ~< W >, which is in an 
explicit contradiction with the direct definition (31). 

Keeping in mind a low-temperature dilute Bose gas, satisfying the Bogolubov theory 
[2-5], we have 

n k = ^coth(^)- I. (32) 



2e k \2TJ 2 

Here the notation 



uj k = ^— + mc 2 
2m 



is used; the Bogolubov spectrum is 



£k — 



2 



and the sound velocity is 



P ^ V^pas 
— st>o — 



m m 

The kinetic energy per particle (31) can be represented as 

K= -p J ^ nk j2nf (33) 

For the local interactions $(r) = $o ( ^( r ) ! the integrand in K behaves as k 4 n k oc $o & t large 
k — > oo, hence, the above integral linearly diverges. However, this by no means implies the 



9 



failure of the Bogolubov theory, but simply reflects a defect in the dilute-gas approximation. 
The Bogolubov theory [2-5] has been formulated for arbitrary interaction potentials $(r), 
but not solely for the local potentials. For a realistic interaction potential $(r), its Fourier 
transform steeply diminishes with increasing k — > oo. Therefore, the integrand in K 
behaves as k 4 rik oc <I>f , which is a rapidly diminishing function of k — > oo. Thus, for realistic 
interaction potentials, no infrared divergence arises in calculating K. Instead of resorting 
to realistic interaction potentials, which would essentially complicate calculations, one can 
model the fast decay of $ fc at large k by the step function 

$ fc = %&(k - k) , (34) 

with 0(/c) being the unit step function. The boundary wave vector k can be chosen as the 
vector, where the phonon behaviour of the Bogolubov spectrum changes to the single- 
particle type dependence. This definition gives ko = 2mc. 

Our goal here is to evaluate the mean kinetic energy (33) by its order of magnitude, in 
order to demonstrate its nonlinear dependence on a s , which makes it principally different 
from the ground-state energy, whose main term is linear in a s . To this end, we resort to 
the simple form (34), assuming that, by the order of magnitude, the integral (33) can be 
evaluated by cutting off its upper limit with the boundary wave vector. Let us also consider 
the zero-temperature case, when the particle distribution (32) becomes 

n k = —r- — (T = 0) . 

2e k 

Then the integral (33) gives 

K = C ^J^ S , (35) 
m v 

C = — (3y/2 - 4) v 7 ^ = 5.505 . 

As is expected, K oc a^ 2 , making it principally different from E oc a s . 

In this way, the Landau-Lifshitz equation (5) is an adiabatic approximation, valid under 
conditions (l)-(3), together with the assumption dS/dm = 0. Incorrectly using this relation, 
outside its region of applicability, leads to wrong paradoxical results. For instance, condition 
(3) does not hold true for the dilute Bose gas, whose interaction potentials do depend on 
particle mass. 

Generally, a parametric differentiation for the statistical average of an operator is given 
by formula (7). If the operator is a self- adjoint integral of motion, it satisfies the zero- 
temperature limit (9). 

The parametric differentiation of the internal energy is represented by the general ex- 
pressions (10), (12), and (13). The mass differentiation of the internal energy is given by 
formula (17). For the dilute Bose gas, the latter takes the form (24), which at zero tem- 
perature reduces to relation (25). In the case of trapped atoms, the mass variation of the 
internal energy leads to Eq. (29). For the homogeneous Bose gas, the mean kinetic energy 
per particle can be evaluated as in Eq. (35). The dependence on the scattering length of 
the kinetic energy K oc cPj 2 and of the internal energy E oc a s is fundamentally different. 



with the constant 
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